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Abstract 

We have analyzed the 241 bursts for which peak counts (C) max exist in the publicly 
available Burst and Transient Source Experiment (BATSE) catalog. Introducing peak 
counts in 1024 ms as a measure of burst brightness B and the ratio of peak counts in 
64 and 1024 ms as a measure of short timescale variability V, we find a statistically 
significant correlation between the brightness and the short timescale variability of 
7-ray bursts. The bursts which are smoother on short timescales are both faint and 
bright, while the bursts which are variable on short timescales are faint only, suggesting 
the existence of two distinct morphological classes of bursts. 

Subject headings: gamma-rays: bursts 

1 INTRODUCTION 

The Burst and Transient Source Experiment (BATSE) has detected more than 600 
7-ray bursts (Meegan et al. 1993). The BATSE catalogue (Fishman et al. 1993) is the 
largest homogeneous sample of 7-ray bursts ever assembled and, as such, constitutes a unique 
resource for studying burst properties. To date, however, no clear correlations between burst 
properties or distinct morphological classes of bursts have been found. 

The Hertzsprung-Russell (H-R) diagram, in which the brightness of a star is plotted 
versus color (the ratio of the flux in two different energy intervals), proved particularly 
useful in classifying stars. A striking feature of 7-ray bursts is the diversity of their time 
histories (see, e.g. Fishman et al. 1992). Guided by the success of the H-R diagram for 
stars and the temporal diversity of 7-ray bursts, we consider an analogous diagram for 7-ray 
bursts in which we plot a measure of burst brightness versus a variability "color" (the ratio 
of the peak flux or the peak counts in two different time intervals) . 



We choose B = (C ,1024 ) max , the expected peak counts in 1024 ms, as our measure of burst 
brightness because 1024 ms is the longest time interval for which peak counts are given in the 
publicly available BATSE burst catalogue, and therefore has the best chance of smoothing 
out effects due to short timescale variability. We choose V = ( (7 64 ) max /(C 1024 ) max , the ratio 
of expected peak counts in 64 and 1024 ms, as our measure of short timescale variability 
because 64 ms is the shortest time interval, and 64 and 1024 ms are the most disparate time 
intervals, for which peak counts are given in the same catalogue. There are many other 
possible measures of variability, and better measures may be found from detailed analysis of 
burst time histories and/or from the acquisition of more data. 

Our measure of variability measures an extreme property of the burst (short timescale 
variability) at a particular moment during the burst (the peak of the burst). Thus it is 
complementary to global measures of burst variability, such as power spectra (Kouveliotou 
et al. 1992) and wavelet analyses (Norris et al. 1992, 1993). 

2 ANALYSIS 

The solid lines in Figure 1 show the time histories of two simulated 7-ray bursts, defined 
by (7 64 , the expected counts in 64 ms. Both bursts have durations tdur = 512 ms; the total 
period shown in each time history, including background, is 1024 ms. For C^ 24 , the expected 
background counts in 1024 ms, we take 2500 counts; this value is typical for the BATSE 
detector (Meegan et al. 1992). For (C' 1024 ) max , the expected peak counts in 1024 ms for the 
simulated 7-ray bursts, we choose 500 counts; this value is typical of the bursts observed by 
BATSE. Since t dur < 1024 ms for both bursts in Figure 1, B = (C 1024 ) max corresponds to 
the total number of counts in the burst, and is therefore a rough measure of burst fluence; 
for bursts with t^nr > 1024 ms, this is no longer true. The time history of the first simulated 
burst has a spike of duration t s J^ e = 64 ms, whereas that of the second is flat. The expected 
peak counts in 64 ms of the two bursts are (C 64 ) max = 200 and 64 counts, respectively. 

We desire B and V as our measures of burst brightness and short timescale variabil- 
ity However, BATSE observes 5 obs = (C 1024 ) max and V ohs = (C 64 ) max /(C 1024 ) max , where 
(C 64 )max and (C 1024 )max are the peak counts in 64 and 1024 ms, respectively. The latter differ 
from the former because the time history observed by BATSE differs from the expected time 
history of the burst due to Poisson fluctuations. The dashed lines in Figure 1 are examples 
of the time histories that BATSE might actually observe for the two bursts. 

Figure 2 (top panel) shows the distribution of burst brightness B ^, s versus short time- 
scale variability V b s for the 201 bursts in the publicly available BATSE catalogue for which 
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both (C 64 ) max and (C 1024 ) max exist. 1 The V" bs-distribution is not independent of 5 Q bs- In 
particular, there are almost no bursts in the upper right quadrant and in a small, triangular- 
shaped region in the lower left hand corner of the diagram. 

By definition, the (C) max protocol selects the largest C 64 and C 1024 during the burst, 
irrespective of the other. Thus (C 64 ) max need not come from the time interval corresponding 
to (C 1024 ) max (nevertheless, one can show that 1/16 < V" b s < 1)- Consider for illustrative 
purposes the limiting case of a flat burst; i.e., a burst for which (C7 64 ) max = C 64 is constant. 
If ^dur < 64 ms, there is only one C 64 during the burst. Although (C 64 ) max = C 64 may differ 
from (C 64 ) max due to Poisson fluctuations, it equals it, on average. If td ur > 64 ms, there are 
several C 64 during the burst, all of which differ from (C) max due to Poisson fluctuations. Since 
the (C) max protocol selects the largest of the C 64 , it is likely that (C 64 ) max exceeds (C 64 ) max 
(see Figure 1). Thus (C 64 ) max is a biased estimator of (C m ) max if ^dur > 64 ms (Meegan 
1993, private communication). Similarly, (C 1024 ) max is a biased estimator of (C 1024: ) maiX . We 
must remove this "Meegan bias" before evaluating the significance of the pattern in the 
(5 obs , V obs )-diagram. 

We have investigated the Meegan bias in detail (Lamb, Boorstein, and Graziani 1993). 
The bias maps the (B, l/)-diagram onto the (.Bobs, Vob s )-diagram. However, the map is not 
invertible because of the Poisson fluctuations of (C 64 ) ma x and (C 1024 ) max away from (C 64 ) maiX 
and (C' 1024 )max, respectively; rather, it tends to shift B and V to larger values of i? b s and 
Vobs, and smears them out. Our studies show that the bias in V bs is zero for a burst whose 
duration tdur < 64 ms, or a burst whose duration tdur is arbitrarily long but which has a 
spike whose duration t s ^ e < 64 ms. The bias in V b s increases until td ur (or t^ c ) reaches 
1024 ms. For longer bursts, the bias increases very slowly because the ratio of the number 
of (C 64 )max samples to the number of (C 1024 ) max samples remains 16. Thus, the bias for 
tdur or t d u^ e = 1024 ms closely approximates the maximum bias for any td ur or fj^°. 

We remove the Meegan bias using an approximate method based on a number of simpli- 
fying assumptions. First, we neglect the bias in B ohs = (C 1024 ) max , which is relatively small. 
Second, we assume that the time history of the burst or of the spike during the burst, is flat; 
in this limiting case, (C m ) m&x = C7 64 and t dur or t d ^ r kc = 641/^ ms. The bias is maximal for 
this model, as noted earlier; any difference that remains between the ^-distributions for faint 
and bright bursts therefore cannot be due to bias. Third, we adopt the bias for the minimum 
of tdur or t d ^ r ke = 64V™ 1 ms and 1024 ms, which closely approximates the maximum bias for 



We have omitted from our analysis four additional bursts (trigger ID numbers 414, 486, 508, and 1346) for 
which the (C 64 ) max or (C 1024 ) max listed in the publicly available BATSE catalogue are erroneous (Howard 
1993, private communication). 
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any tdur or f^ . Fourth, we assume that the bias maps a point in the (B, V)-diagram onto 
a unique point in the (-B obs , V bs)-diagram; i.e., that the bias map is a 5-function. This is 
only approximately true because of the variances of (C 64 ) max and (C 1024 ) max about (C 64 ) max 
and (C 1024 ) max . These four simplifying assumptions enable us to invert the bias, and map 
bursts in the (.Bobs, V" b s )-diagram onto points in the (B, F)-diagram. 

Figure 2 (middle panel) shows contours of constant V in the (B b s , K>bs) -diagram, taking 
for the expected background counts in 1024 ms C^ 2A = 2500 counts, a value which is typical 
for the BATSE detectors. The bias is zero along the top of the diagram. Elsewhere, the 
bias is least in the upper right hand corner and greatest in the lower left hand corner of 
the diagram. This pattern arises partly from the time variability of the bursts, which is 
greatest along the top of the diagram and least along the bottom, and partly from the 
relative variances in the total counts + (C 64 )max and C 1 ^ 24 + (C ,1024 ) max , which are least 
in the upper right hand corner of the diagram and greatest in the lower left hand corner. 

Bursts that lie below the V = 1/16 contour (shown as a solid line) in the lower left hand 
corner of the (B obs , Vo bs )-diagram would map onto points below V — 1/16 in the (B, V)- 
diagram, which is unphysical. This is due to our assumption that the burst time history is 
flat, for which the bias is maximal, and to our 5-function approximation for the bias map. 
We therefore place a "floor" at V — 1/16 in the (B, VQ-diagram, and do not allow points 
mapped from the (B ohs , V obs )- diagram to go below it. 

3 RESULTS 

Figure 2 (bottom panel) shows the resulting distribution of bursts in the (B, F)-diagram. 
The statistical errors in B and V are a function only of location in the (B, VQ-diagram; we 
show them for six representative locations. The empty triangular region at the lower left 
hand corner of the (B b s , Vob s )-diagram has disappeared; it may therefore be the result of 
Meegan bias. In contrast, the presence of bursts in the upper left hand quadrant and the 
absence of bursts in the upper right hand quadrant of the diagram remain. This means that 
there is a lack of bursts that are bright (in 1024 ms) and have a bright, short spike (either 
during the bright 1024 ms or elsewhere during the burst). 

To evaluate the significance of the pattern, we use three \ 2 tests for binned data which 
address the question of whether or not two distributions have similar shapes, irrespective 
of scale, i.e., irrespective of the total number of objects in each histogram. The first is 
the "T-test" of Eadie et al. (1971), in which the scales of the two histograms are set at 
their maximum likelihood best-fit values. The second is the two-histogram test described 
in Numerical Recipes (Press et al. 1986), with the histogram having the smaller number of 
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TABLE 1 



Histograms of faint, Bright, 
Smooth, and variable Bursts 



Two Bin Three Bin 

Sample Bin 1 Bin 2 Bin 1 Bin2 Bin 3 



Faint 


135 


27 


126 


19 


17 


Bright 


38 


1 


38 





1 


Smooth 


131 


33 


81 


73 


10 


Variable 


36 


1 


28 


8 


1 



objects scaled up so that both histograms have the same number of objects. The third is an 
adaptation of the previous test in which the scale of the histogram with the smaller number 
of objects is allowed to vary so as to minimize x 2 

We divide the bursts into those which are faint (log B < 3.28) and bright (log B > 3.28), 
and construct histograms having two and three equal logarithmic bins in variability between 
logl^ = —1.2 and 0. We also divide the bursts into those which are smooth (log V < —0.8) 
and variable (log V > —0.8) on short timescales, and construct histograms having two and 
three equal logarithmic bins in brightness between log B = 2 and 4.7. Table 1 gives the 
numbers of bursts in each bin of these histograms, and Table 2 lists the results of the three 
tests. 2 The Q-values range from 2.2 x 10~ 2 to 3.5 x 10~ 6 , comparing the variability of faint 
and bright bursts, and from 1.5 x 10~ 2 to 7.8 x 10" 5 , comparing the brightness of smooth 
and variable bursts. 

The disparity between the results of the three tests is due primarily to the small number 
of bursts in some histogram bins (see footnote 2). We therefore also evaluate the significance 
of the correlation in the (B, l / )-diagram using a statistical test which compares the mean 
variabilities of faint and bright bursts, or the mean brightnesses of smooth and variable 
bursts, and evaluates the significance of the differences between the means using a t-test 
(Press et al. 1986). This test is particularly suitable to the situation at hand because it 
implicitly takes into account the uncertainty in the location of each burst in the (B, V)- 
diagram and, unlike x 2 tests, does not require binning the data, which can necessitate the 



2 The results are similar for histograms with more bins. However, the results are increasingly suspect because 
the number of bursts becomes small in many histogram bins, necessitating the use of Poisson rather than 
Gaussian statistics and invalidating the use of the x 2 distribution for calculating statistical significance. 
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TABLE 2 



Comparison Between faint and Bright, and Smooth and variable Bursts 







Faint vs. 


Bright 




Smooth vs. Variable 


Statistical 




2 Bins 


3 Bins 




2 Bins 3 Bins 


Test 


x 2 


Q-value 


X 2 Q-value 


x 2 


Q-value \ 2 Q-value 


T-Test a 


5.2 


2.2 x 10~ 2 


8.3 1.6 x 10"2 


6.5 


1.0 x 10~ 2 8.4 1.5 x 1(T 2 


Fixed Scale b 


12 


4.2 x 10" 4 


25 3.5 x 1(T 6 


16 


4.9 x 10" 5 10 6.4 x 10~ 3 


Variable Scale c 


12 


5.1 x 10~ 4 


25 3.7 x 10" 6 


16 


7.8 x 10~ 5 8.9 1.2 x 10~ 2 



a Two-histogram "T-test" (Eadie et al. 1971). 
b Two- histogram test (Press et al. 1986). 

c Adaptation of two-histogram test in Press et al. (1986) (see text). 



introduction of Poisson statistics when the number of bursts is small. 

We find that the faint (B < 1900) bursts have a mean short timescale variability log V = 
-0.97 ± 0.026 whereas the bright (B > 1900) bursts have TogV = -1.05 ± 0.028. Thus the 
difference of the means divided by their combined variance is 2.2, corresponding to a Q- value 
of 3.3 x 10~ 2 . We find that the smooth bursts have a mean brightness log 5 = 2.99 ± 0.036 
whereas the variable bursts have log B = 2.67 ± 0.059. Thus the difference of the means 
divided by their combined variance is 4.5, corresponding to a Q-value of 4.6 x 10~ 5 . 

We conclude that the difference in short timescale variability between faint and bright 
7-bursts, and the difference in brightness between bursts that are smooth and variable on 
short timescales are significant. 

We have checked our results in two ways. First, we determined the number of bursts 
in our study which have gaps in their time histories, due to malfunctioning of the Compton 
Observatory tape recorders. Using the comments table in the publicly available BATSE 
catalogue, we find that 25 of the 48 bursts with data gaps are among the bursts we have 
used. Of these 25, 21 are faint and 4 are bright bursts; three of the faint bursts and none 
of the bright bursts are variable on short timescales. We conclude that gaps in the time 
histories of the bursts cannot account for our results. 

Second, we inspected the time histories of the 30 individual bursts out of the 201 bursts 
in our sample for which time histories are publicly available. 3 Of these, 17 are faint and 13 



3 Two are in Fishman et al. (1992a) and ten are in Fishman et al. (1992b); eighteen others are among the 
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are bright bursts; five of the faint bursts and one of the bright bursts are variable on short 
timescales. In all cases, we were able to verify from inspection of the time history that the 
location of the burst in the (B, V^-diagram is correct. 

4 DISCUSSION 

The correlation between the brightness B and the short timescale variability V of 7-ray 
bursts might imply strong source evolution, or it might reflect the existence of two distinct 
classes of 7-ray bursts. One class might be bursts that exhibit a range of short timescale 
variability (—1.2 < \ogV < 0) and are faint, and the other might be bursts that are smooth 
on short timescales (logV < —0.8) and bright. Alternatively, one class might be bursts that 
are smooth on short timescales (log V < —0.8) and range from faint to bright, and the other 
might be bursts that are variable on short timescales (logV > —0.8, corresponding to t dur 
or t s J^° ^ 0.3 s) and faint. Clearly, other decompositions are also possible. 

Using only the correlation between B and V presented here, we cannot distinguish 
between these various possibilities. However, in subsequent papers, we report statistically 
significant evidence that the correlation between B and V reflects the existence of two distinct 
morphological classes of 7-ray bursts: Type I bursts, which are smooth on short timescales 
0.3 s) and range from faint to bright, and Type II bursts, which are variable on short 
timescales and faint. Type I bursts also have longer durations and softer spectra (Lamb 
and Graziani 1993a), and a flatter brightness distribution (Lamb and Graziani 1993b) than 
do Type II bursts. The dashed lines in the top and bottom panels of Figure 2 correspond 
to \ogV = —0.8, the cut in variability which separates the two classes. Because we have 
removed the Meegan bias assuming that burst time histories are flat, for which the bias is 
maximal, it is likely that we have placed some bursts which are variable on short timescales 
in the smooth class, but not vice-versa. 
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here. We thank Tom Loredo for many valuable discussions about statistical methodology, 
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of expected counts. We thank Josh Boorstein for help in investigating the nature of the bias 
and Peter Freeman for help in obtaining burst time histories. CG acknowledges the support 
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Figure 1: Time histories of two simulated 7-ray bursts, each lasting tdur = 512 ms. (left 
panel) Burst with a spike in (7 64 lasting t s J^° = 64 ms. (right panel) Burst with C u equal 
to a constant. The solid histograms show (7 64 , the expected counts per 64 ms; the dashed 
histograms show C 64 , the observed counts per 64 ms. Also labeled are (C' 64 ) max , the expected 
peak counts per 64 ms, and (C 64 ) max , the observed peak counts per 64 ms. 

Figure 2: (top panel) Distribution of 201 bursts in the (5 bs, V obs )-diagram. (middle panel) 
Contours of constant V in the (i? bs, K>bs) _ diagram. (bottom panel) Distribution of 201 
bursts in the (B, F)-diagram. The filled and open circles in the (-Bobs 5 V b s )- and (B, 
V)-diagrams denote bursts detected by the 1024 ms trigger and by the 256 or 64 ms trig- 
gers, respectively. The solid lines in the (B^, V b s )- and (B, l/)-diagrams correspond 
to V — 1/16, the minimum possible variability, while the dashed lines correspond to 
log V = —0.8, the cut in variability which separates the bursts into two distinct classes. 



This manuscript was prepared with the AAS WGAS IATgX macros v2.2 



9 



